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We calculate the linear and non-linear susceptibilities of periodic longitudinal chains of hydrogen 
dimers with different bond-length alternations using a diffusion quantum Monte Carlo approach. 
These quantities are derived from the changes in electronic polarization as a function of applied 
finite electric field - an approach we recently introduced and made possible by the use of a Berry- 
phase, many-body electric-enthalpy functional. Calculated susceptibilities and hyper-susceptibilities 
are found to be in excellent agreement with the best estimates available from quantum chemistry 
- usually extrapolations to the infinite-chain limit of calculations for chains of finite length. It is 
found that while exchange effects dominate the proper description of the susceptibilities, second 
hyper-susceptibilities are greatly affected by electronic correlations. We also assess how different 
approximations to the nodal surface of the many-body wavefunction affect the accuracy of the 
calculated susceptibilities. 

PACS numbers: 36.20.-r, 77.22.-d, 71.15.-m, 31.15.A- 



I. INTRODUCTION 

The linear and non-linear longitudinal dielectric sus- 
ceptibilities of periodic linear chains of hydrogen dimers 
have been in recent years the subject of many density- 
functional theory and quantum chemistry studies 
li,.5, 6, 7,.8, Ji,[ia,[lll,[l2, 13]. When dielectric properties 
are considered, these model systems present appealing 
similarities with real polymeric chains. In particular, in 
both systems the second hyper-susceptibilities are very 
large. Despite the apparent simplicity of such systems, 
obtaining reliable results for susceptibilities and hyper- 
susceptibilities requires particular care [U- Density- 
functional theory approaches using the local density 
(LDA) or generalized-gradient (GGA) approximations 
strongly overestimate dielectric response, up to several 
orders of magnitude in the case of the second hyper- 
susceptibility, and more sophisticated approaches need 
to be considered [isj . Correlated quantum chemistry 
approaches, such as M0ller-Plesset (MP) and coupled- 
cluster (CC), provide benchmark results, although care- 
fully converged results often require the most sophisti- 
cated levels of approximations (MP4,CCSD(T)) and ex- 
tended basis setsQ]- 

Recently, we have introduced a novel scheme for the 
treatment of finite static homogeneous electric fields in 
periodic quantum Monte Carlo (QMC) calculations [l6|. 
The use of periodic boundary conditions (PBCs) permits 
to avoid the extrapolations to the infinite chain limit, 
which would be even more cumbersome in this case due 
to the presence of QMC statistical errors. When applied 
to the hydrogen chain model, this approach provides lin- 
ear susceptibilities in excellent agreement with the results 
from the most accurate quantum chemistry calculations. 
Since it has been shown that, in contrast to the linear re- 
sponse case, electronic correlations account for a large ex- 
tent to second hyper-susceptibilities 01, it is of particular 



interest to assess the accuracy and reliability of QMC for 
the calculation of these higher-order responses. Indeed in 
QMC, in the limit of a correct nodal surface approxima- 
tion, electronic correlations are treated exactly and basis 
sets approximations are easily controlled. Moreover, due 
to the favorable scaling of QMC approaches with respect 
to the number of electrons (cubic, or even linear 17]), 
this approach could represent a very viable strategy for 
accurate, predictive calculations of non-linear dielectric 
properties in realistic systems. 

Here, we calculate with diffusion quantum Monte 
Carlo (DMC) the linear polarizability and second hyper- 
polar iz ability per II2 unit (i.e. the susceptibilities) for the 
paradigmatic case of periodic linear chains of hydrogen 
dimers, where the dimers have a fixed length of 2.0 a. u. 
and we consider three different bond-length alternations 
of 2.5 a. u. 3.0 a. u. and 4.0 a. u. respectively. 

The paper is organized as follows: Section |ll] briefly 
describes our methodology to treat finite electric fields in 
periodic QMC calculations. In Section ITTTl we summarize 
the technical details of the calculations and investigate 
the degree of convergence with respect to the different 
approximations used. The final results are reported and 
discussed in Section ITVl while in Section|V]we investigate 
the errors induced by different choices for the many-body 
nodal surface. Conclusions, and perspectives for future 
work are drawn in Section [VII 



II. METHOD 

We briefly summarize here the approach we first intro- 
duced in Ref. [16^ . We consider a system of N interacting 
electrons in a periodic cell of size L (to simplify the nota- 
tion we describe here the one-dimensional case, but the 
extension to higher dimensions is straightforward). The 
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normalized many-body wavefunction obeys PBCs 
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for every i. On the other hand, the position operator 
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does not satisfy PBCs and is unbounded from below; thus 
it cannot be used (and would be ill-defined) to calcu- 
late expectation values, say, of the potential generated 
by a constant electric field. On the other hand, the 
modern theory of the polarization[l^[l^, [l^l (MTP) and 
its many-body generalization [IH provide definitions for 
the relevant observables which remains compatible with 
PBCs. Indeed, the polarization P of a many-body sys- 
tem can be obtained from the single-point Berry phase 
(21I formulation 
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where ft is the size of the cell (in Id is equal to L), e is the 
elementary charge, G — 2tt/L, and z will be termed here 
as the complex polarization. The definition of Eq. ([3|) 
coincides, in the thermodynamic limit, with the exact 
many-body observable [21|, but it remains well-defined 
for any finite L. Although in Eq. ([3|) the polarization is 
defined modulo —eL/Q, the quantities which are experi- 
mentally observed correspond to differences of polariza- 
tion and are well defined within the MTP. 

Having a valid definition for the many-body polariza- 
tion, we can operate a Legendre transform of our total 
energy functional between the conjugate variables of elec- 
tric field £ and polarization P, and obtain the ground- 
state wavefunction of a periodic, extended system in the 
presence of an electric field £ from the minimum of the 
generalized electric enthalpy [H [13, HI [H HI] 



(5) 



where [^] is the energy functional for the unperturbed 
Hamiltonian iJ". The direct minimization of an electric 
enthalpy functional has been firstly introduced for tight- 
binding Hamiltonians in Ref. [2J] and in the context 
of density functional theory in Refs. [13, [13|- On the 
other hand, DMC and variational Monte Carlo (VMC) 
(with variance minimization approaches) require the use 
of Hamiltonian operators which are both local and Her- 
mitian. This requirement can be satisfied by noting that 
a local Hermitian operator can be identified from the 
minimum condition for the electric enthalpy 



6F 

m 



(6) 



with A the appropriate Lagrange multiplier. This strat- 
egy has already been used for the case of single-particle 



Slater determinants and model Hubbard Hamiltoni- 
ans [1^. This local Hermitian operator shares the same 
ground-state wavefunction of the electric enthalpy func- 
tional, and the ^' that minimizes Eq. ([5]) is also the 
ground-state for the many-body Hamiltonian 



H{z) = H" 
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It should be noted that Eq. ([7]) defines a self-consistent 
many-body Hamiltonian, since it's an operator that de- 
pends on the complex polarization 2, and thus on the 
expectation value of through Eq. Q . Self-consistency 
is indeed the price to be paid to recast the variational 
principles of the electric-enthalpy functional into an op- 
erator equation. 

Due to the self-consistent nature of the operator H{z) 
defined in Eq. ([7]), the ground state in the presence of an 
electric field must be found through an iterative proce- 
dure. We start from a first value zi for z, e.g. as found in 
the single-particle calculations or in the many-body trial 
wave function $t; the local Hamiltonian H{zi) is then 
constructed. DMC evolution using this operator leads to 
a new expectation value for z, called Z2, which in turn 
determines a second Hamiltonian H{z2). In the absence 
of stochastic noise, this process could be iterated to con- 
vergence: 



Z2 



Z3 



to find the fixed point of the complex-plane map 

f{Zi) = Zi+l. 



(8) 



(9) 



Since the Monte Carlo procedure introduces a statisti- 
cal error in every estimate of z.^, the map / becomes a 
stochastic function in the complex plane. If / can be ap- 
proximated close to its fixed point as a linear function in 
the complex plane, it has been shown[l6] that its aver- 
age over a sequence of {z^} provides the best estimate for 
the fixed point itself. The validity of this approximation 
must be checked numerically. Once the fixed point z is 
obtained, the corresponding polarization is then obtained 
via Eq. Q. 



III. COMPUTATIONAL DETAILS 

We investigate in this study dielectric properties of pe- 
riodic linear chains of hydrogen dimers along the longi- 
tudinal direction. The chosen dimer bond length is 2.0 
a. u. and bond-length alternations of 2.5, 3.0 and 4.0 a. 
u. are considered. Because of symmetry, these systems 
have vanishing first hyper-susceptibility. We treat these 
systems through periodic boundary conditions. We first 
perform HF calculations using an orthonormal unit cell 
with an edge size of 20 a. u. in the directions perpen- 
dicular to the chain and containing one single H2 unit. 
We sample the one-dimensional Brillouin zone along the 
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chain direction with a set of Nk equally spaced k-points. 
Then, for QMC calculations the corresponding supercell 
containing Nk hydrogen dimers, sampled at the F-point, 
is considered. The HF calculations are performed us- 
ing the PWSCF package from the Quantum-ESPRESSO 
distribution (28|. The wave- functions are expanded in a 
plane-wavebasis set with a cutoff of 50 Ry. Finite electric 
fields are implemented through the algorithm described 
in Ref. jl^l ■ To describe the Coulomb potential of the H 
nucleus we use a norm-conserving pseudopotential. For 
this same pseudopotential we already showed in Ref. [l6l | 
that the calculated DMC polarizability of the isolated H 
atom is in perfect agreement with the theoretical value. 

The convergence of the Berry-phase polarization to the 
thermodynamic limit is of the order with respect to 
the supercell dimension L; this can be understood on 
purely geometrical grounds (see Ref. 29]). Thus, we de- 
termine here the k-point sampling Nk for the primitive 
unit cell, and equivalently the length of a supercell sam- 
pled at the F point, for which the sampling/geometrical 
errors in the calculated linear susceptibility and second 
hyper-susceptibility become smaller or negligible with re- 
spect to the magnitude of the statistical errors due to the 
DMC procedure. We investigate this convergence within 
the HF approach, focusing on chains with a bond-length 
alternations of 3 a. u. , for which also extrapolations to 
the infinite chain limit for both the polarizability and 
hyper-polarizability have been reported for several basis 
sets in Ref. 0. We consider two values for Nk, 10 and 20, 
and use two applied electric fields of magnitude 0.003 and 
0.02 a. u. respectively. We then fit the dipole moments 
fi per H2 unit with 

=aS + l^S^, (10) 
6 

where £ is the electric field and a and 7 are the linear 
susceptibility and the second hyper-susceptibility per H2 
unit. We report in Tab. [T] the values for a and 7 we 
calculated using PBCs, together with the extrapolations 
for the most complete basis sets presented in Refs. [ll[3> 
showing very good agreement with these published cal- 
culations. As we want to focus here on the non-linear 
susceptibility we will consider in the following supercells 
of 10 dimers each, i.e. corresponding to Nk = 10 sam- 
pling. 

The calculated HF wavefunctions are then interpolated 
with splines and imported into the CASINO VMC and 
DMC codes [30l|. In VMC the iV-body wavefunction 
^VMC is defined as 

*VMc(ri, •■,rAr) = exp{J{ri, .., rAr))— 1= det (V'i..'0jv) , 

(11) 

where {r} are the electron positions, {i/j} are the HF 
single-electron wavefunctions and J is the Jastrow fac- 
tor. In this work we use the Jastrow factor formulation 
introduced in Ref. [3l| . which depends only on the dis- 
tances |rj — r^l and jr^ — R/|, where {R} are the ionic 
positions. We use an expansion to the sixth order for 



both spin components, and two different sets of parame- 
ters for inequivalent H ions. Indeed, when a longitudinal 
field is applied, the two H atoms in the unit cell are no 
longer equivalent. To optimize the Jastrow parameters 
we apply a variance minimization scheme [s^. For VMC 
simulations we adjust the time-step in order to assure an 
acceptance ratio of ~50%. It should be noted that we 
use a spherical symmetrical Jastrow factor that depends 
only on the relative distance between electrons and be- 
tween electrons and nuclei. Such form does not allow, in 
the VMC procedure, for significant changes of the expec- 
tation value of z (defined in Eq. ([3])) with respect to its 
original HF value. Moreover, the variance minimization 
scheme used here is less sensitive to lon g ra nge properties 
than energy minimization schemes [s^ |33|| 

In the following, we consider a hydrogen chain with 
a bond- length alternation of 2.5 a. u. , as done in Ref. 
[l^ . We report in Tab.HIlthe estimates for the total elec- 
tric enthalpy and z complex polarization obtained in the 
presence of an electric field of 0.003 a. u. and for several 
consecutive steps of the Jastrow factor optimization. As 
expected, we note a significant lowering of the total elec- 
tric enthalpy due to the term in Eq. while the 
polarization remains almost unchanged. 

Once the Jastrow factor has been optimized, the result- 
ing A^-body wavefunction ^vMC is used as trial wave- 
function for the DMC simulations, which are based on 
importance sampling. In all the DMC calculations we 
use a time step of 0.02 a. u. assuring an acceptance ratio 
greater than 99.7%. Indeed, in the DMC simulations the 
walkers are distributed as \E'vmc^', where is the cor- 
rect A^-body ground-state wavefunction compatible with 
the nodal surface defined by ^'vmc- The expectation 
value O of an operator O which does not commute with 
the Hamiltonian and depends only on the electron posi- 
tions |r| is obtained through a forward-walking proce- 
dure 

where Nt is the number of walkers at time-step r and 
Tj^r-At corresponds to the ancestor configuration of the 
walker j at an interval At back in imaginary time. To 
determine the projection time At we start from a 5'vmc 
calculated for an electric field of 0.003 a. u. and switch 
off the electric field £ in the following DMC simulation. 
Then we calculate the expectation values of z for increas- 
ing values for At and through Eq. ^ the corresponding 
dipole moments. These are expected to decay exponen- 
tially toward a small value. Indeed constraining 5* to the 
nodal surface of which is calculated for a non- 

zero electric field, prevents the complete relaxation of the 
dipole moment when the field is switched off in the DMC 
simulation. In Tab. lIIII we display the decay of the dipole 
moment for At ranging from to 1000 time steps. In this 
study we then use a At corresponding to 1000 time steps 
which assures the convergence of the estimated values of 
z with respect to the statistical errors of our DMC simu- 
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lations. We do not report in the calculated DMC errors 
the residual bias due to the use of the finite projection 
time A. 

Last, to verify the linearity of the map / in Eq. ([9|) 
in the vicinity of its fixed point we perform two differ- 
ent sequences of DMC simulations, as described in Eq. 
([8]), starting with the z value given by HF in an applied 
electric field of 0.003 a. u. . While both series have a 
population of 2560 walkers, the first one is composed of 
20 DMC simulations of 40000 time steps each, while the 
second one is composed of 6 DMC simulations of 120000 
time steps each. Therefore, larger error bars for each es- 
timate z are found in the first series. We expect that if 
/ is linear the calculated average values J. should be the 
same for the first and second series, within the statisti- 
cal error. Indeed, we report in Tab. IIVI the average z 
together with the corresponding dipole moments found 
for the two series. Both quantities are found to be equal 
within the statistical errors. Thus, in the following DMC 
simulations we use runs of 40000 time steps each and 
with a population of 2560 walkers. 



IV. RESULTS AND DISCUSSION 

Using the computational scheme and optimal parame- 
ters described in the previous sections we consider linear 
hydrogen chains with bond-length alternations of 2.5, 3.0 
and 4.0 a. u. , intermolecular distances of 2.0 a. u. de- 
scribed with 10-dimer supercells in PBCs, to calculate 
the DMC best estimates for the fixed-point z at two dif- 
ferent values of the applied electric field. We report in 
Tab. |V] these estimates z for the complex polarization, 
together with the corresponding dipole moments per H2 
unit. Since these systems are centrosymmetric, we con- 
sider the dipole moment to be zero in the absence an 
electric field. For larger bond-length alternations it is 
possible to access larger field intensities without runaway 
solutions in the HF calculations [22]. 

Then, we fit the calculated moments with the expan- 
sion of Eq. (llOp . and obtain the DMC estimates for the 
linear susceptibility a and second hyper-susceptibility 
7. The linear susceptibilities are reported in Tab. IVII 
together with the quantum chemistry extrapolations of 
Ref. We determine a with statistical errors vary- 

ing from 1.5 % to 0.4 % for the 3.0 a. u. and 4.0 a. u. 
bond-length alternations, respectively. The higher preci- 
sion in the latter case is due to the application of higher 
electric fields. The reported M0ller-Plesset susceptibili- 
ties calculated for the levels of approximation MP3 and 
MP4 and for the (6)-31G(*)* and (6)-311G(*)* basis sets 
show a good degree of convergence with the level of ap- 
proximation (from MP3 to MP4) and a slightly worse 
degree of convergence when using a more complete basis 
set ( from (6)-31G(*)* to (6)-311G(*)*). Similar results 
to the MP4 case are obtained using the coupled cluster 
CCSD(T) method, as illustrated in Ref. [H]. Our DMC 
estimates for a are in excellent agreement with the MP4 



values calculated for the basis set (6)-311G(*)*. In Ref. 
[T6| . the same good agreement was found for the chain 
with a bond-length alternation of 2.5 a. u. and using a 
trial wavefunction obtained from a density- functional cal- 
culation. 

We show in Tab. IVIII the calculated DMC estimates 
for the second hyper-susceptibility 7, together with the 
quantum chemistry extrapolations of Ref. [2]. We de- 
termine 7 with statistical errors varying from 6.8 % to 
3.6 % for the 3.0 a. u. and 4.0 a. u. bond-length alterna- 
tions, respectively. Values from MP quantum chemistry 
extrapolations are available only for the 3.0 a. u. bond- 
length alternation case. These show slow convergence in 
the MP series (MP3 to MP4) and are still not converged 
with the basis set (from (6)-31G(*)* to (6)-311G(*)*). 
However, in Ref. [2] it has been shown that 7 calculated 
through MP3 and MP4 increases while the basis set is be- 
coming more complete. Therefore, the best MP4 value 
should be taken as a lower limit. In fact, the 7 coefficient 
calculated through DMC is found to be higher than the 
MP4 counterpart by a factor of ~20%. 

By addressing the difference between the susceptibili- 
ties calculated with HF and those calculated with DMC, 
we can investigate the relevance of electronic correla- 
tions. We report in Tab. IVIIII these differences with 
respect to the corresponding DMC susceptibilities and 
hyper-susceptibilities. Electronic correlations lower the 
value of linear susceptibilities a and increase that of sec- 
ond hyper-susceptibilities 7. As already shown in Ref. 
0, the effects of correlations are large for the second 
hyper-susceptibility 7. Indeed, they account for almost 
half of the calculated DMC value for the chain with 
a bond-length alternation of 2.5 a. u. . As the bond- 
length alternation increases the correlation contribution 
becomes smaller but still remains considerable. There- 
fore, when calculating non-linear susceptibilities in such 
systems, an adequate treatment of electronic correlations 
is mandatory and QMC approaches are particularly ap- 
pealing since in the limit of the exact nodal surface they 
treat correlations exactly. 

It is worth investigating the localization of the iV-body 
wavefunction ^E* along the longitudinal chain direction. 
The MTP provides a definition of the localization spread 
which remain valid also within PBCs[33|- This de- 
pends solely on the complex number z which we esti- 
mated in DMC, and the following holds 
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We report in Tab. llXl the DMC spreads cr^ for the bond- 
length alternations and electric fields addressed. We note 
that as the bond-length alternation increases the sys- 
tem becomes more localized. Indeed in the limit of large 
bond-length alternations, the system becomes composed 
of isolated hydrogen dimers. The opposite behavior is 
observed when the electric field is increased. Indeed an 
applied electric field closes the electronic gap, and this 
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determines an increase of the electronic spread as shown 
in Ref. [H. 

Finally, we address the converge of the calculated DMC 
linear susceptibilities and second-hyper susceptibilities 
with respect to the size of the supercell used in the sim- 
ulation. We consider a bond length alternation of 4 a. u. 
and we perform an additional calculation for a supercell 
consisting of 20 H2 units. We use the same calculation 
parameters of the previous DMC simulations. We see 
from the figures reported in Tab. |X] that, in the limit 
of the statistical error, many-body effects play a minor 
role. However, a full assessment of the convergence of the 
second-hyper susceptibility with respect to the supercell 
size would require significant smaller error-bars. 

V. NODAL SURFACE APPROXIMATION 

In this section we address the dependence of the cal- 
culated DMC dipolc moments with respect to the nodal 
surface through the choice of the trial wavefunction. To 
estimate the magnitude of the error due to an approxi- 
mate nodal surface, we consider the hydrogen chain with 
a bond- length alternation of 2.5 a. u. and proceed in this 
way: first we calculate HF trial wave-functions applying 
electric field intensities of 0.0, 0.003 and 0.01 a. u. then, 
the following VMC and DMC simulations are performed 
for an electric field of 0.003 a. u. and the fixed points 
z and dipole moments are compared. The discrepan- 
cies between the DMC values calculated starting from 
HF wavefunctions obtained consistently for an electric 
field intensity of 0.003 a. u. and those from wavefunc- 
tions obtained at fields of 0. and 0.01 a. u. are due to 
the different nodal surfaces. These values are displayed 
in Tab. IXII together with the corresponding dipole mo- 
ments. We note that for our system the imaginary part 
of z is most affected by the choice of the starting trial 
wavefunction. As expected, the dipole moment is found 
to be smaller when starting from an HF wavefunction 
presenting a smaller dipole moment and becomes larger 
when the HF wavefunction has a larger one. However, 
this discrepancy is not very large: when starting from 
a 0. a. u. dipole moment a difference of only ~ 7% is 
observed. This discrepancy increases up to ~ 20% (i.e. 
three times larger, but still much smaller than the differ- 
ence with the initial HF value) when starting from the HF 
wavefunction with a dipole moment of 0.6055 a. u. . We 



can conclude that the error in the DMC dipole moments 
is negligible when using trial wavefunctions calculated for 
the same values of the electric field. The same conclusion 
was drawn in Ref. ^1^, where it was shown that when 
starting from LDA trial wavefunctions, which strongly 
overestimate the linear susceptibilities, the DMC linear 
susceptibility was still in excellent agreement with the 
best quantum chemistry results. This could be due to 
the resemblance with genuine one dimensional systems 
for which the nodal surface is fixed by symmetry. 
VI. CONCLUSIONS AND PERSPECTIVES 

We have shown how it is possible to calculate linear 
and non-linear susceptibilities of periodic systems using 
DMC, obtaining for the cases studied statistical errors 
lower than 1.5% for the linear susceptibilities and 7% 
for the non-linear susceptibilities. The calculated values 
are in excellent agreement (when available) with the re- 
sults obtained from the most accurate quantum chem- 
istry approaches, but do not suffer from basis set er- 
rors or extrapolations to the infinite limit. These results 
and approach can be expected to be used as reference 
for testing novel first-principles approaches for the eval- 
uation of dielectrics properties. Although DMC calcu- 
lations of susceptibilities are computational demanding, 
they can benefit, in terms of computational speed, from 
the use of order- methods based on localized orthogo- 
nal [13, nil or nonorthogonal basis sets [10, Slj . There- 
fore, it would be of great interest to determine the ac- 
curacy of the calculated dielectric properties when such 
approximations are used. DMC estimates of linear and 
non-linear dielectric properties could be particularly im- 
portant for systems for which less expensive approaches, 
such as density-functional theory, fail. Polymers and con- 
jugated organic molecules are a prototypical class of these 
systems, of great theoretical and practical interest. 
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Nk = 10 Nk= 20 (6)-31G(*)* (6)-311G(*)* 
a 28.5 28.5 28.5 28.6 
7 56.0 57.1 55.1 56.7 

TABLE I: Hartree-Fock linear susceptibility a in a. u. per H2 
unit and second hyper-susceptibility 7 in 10^ a. u. per H2 unit 
for the periodic linear chain of H2 dimers (d=2.Q a. u. ) with 
a bond-length alternation of 3 a. u. , calculated using PBCs 
and meshes of Nk = 10 and A^^ = 20 equally spaced k-points. 
The quantum chemistry extrapolations to the infinite limit, 
obtained with basis sets (6)-31G(*)* and (6)-311G(*)*, are 
taken from Refs. P, Q]. 

#VMC run F 

-10.6855 ± 0.0011 0.6118 ± 0.0004 0.1443 ± 0.0004 TABLE III: DMC dipole moment fi per Ha unit for a lin- 

1 -11.1162 ± 0.0015 0.6128 ± 0.0013 0.1453 ± 0.0013 ear periodic chain (10-dimer supercell) with a bond-length 

2 -11.1608 ± 0.0014 0.6128 ± 0.0014 0.1472 ± 0.0014 alternation of 2.5 a. u. , as a function of the forward-walking 

3 -11.1582 ± 0.0014 0.6141 ± 0.0014 0.1466 ± 0.0014 P''°j*''^*^°'^ ^^'^ ^^"^ ^^'^ ^"^^ wavefunction used has 

been determined with an applied electric field of 0.003 a. u. 



At (time steps) 


H (a. u. ) 





30.8 ± 1.3 


100 


22.8 ± 1.4 


200 


17.6 ± 1.5 


300 


13.9 ± 1.6 


400 


11.2 ± 1.7 


500 


9.2 ± 1.8 


600 


7.8 ± 2.0 


700 


6.8 ± 2.1 


800 


6.1 ± 2.3 


900 


5.6 ± 2.4 


1000 


5.4 ± 2.6 



rnA T>T n TT i-M i • ii, 1 J 1 1 • i' No elcctric ficM Is appllcd durlug thc DMC slmulatious, cach 

lABLE 11: Electric enthalpy r and complex polarization z , 

■u. ■ J ■ A i- Tnvj/-i c • J- u J involving 1000 walkers and 36000 time steps, 

obtained m 4 consecutive VMC runs tor a periodic hydrogen ^ ^ 

linear chain of H2 dimers (ci=2.0 a. u. ) with a bond-length 
alternation of 2.5 a. u. and an applied electric field of 0.003 
a. u. . The supercell contains 10 dimers. In the 0-th run no 
Jastrow term is used and the simulation has a population of 
3.2*10^ walkers. In runs 1-3 the Jastrow factor is optimized 
starting from the parameters used in Ref . [l^ , and the simu- 
lations have a population of 3.2*10^ walkers. 
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# time-steps # runs ?Si{z} ^{z} l-i (a. u. ) 

40000 20 0.6122 ± 0.0007 0.1345 ± 0.0016 0.1549 ± 0.0020 

120000 6 0.6130 ± 0.0008 0.1334 ± 0.0015 0.1545 ± 0.0018 

TABLE IV: Real and imaginary part of the estimated fixed 
point z for the complex polarization and corresponding dipole 
moment ^ per H2 unit, calculated for a linear periodic chain 
(10-dimer supercell) with a bond-length alternation of 2.5 a. 
u. and an applied electric field of 0.003 a. u. , for two iterative 
series with time-steps) time-steps per single run, and 
runs) total number of DMC runs. 



L(a. u. ' 


) £ (a. u. ' 


) ^{z} 


Q{z} 


fi (a. u. ) 


2.5 


0.003 


0.6127 ± 0.0008 


0.1341 ± 0.0013 


0.1544 ± 0.0018 


2.5 


0.01 


0.3736 ± 0.0020 


0.4320 ± 0.0034 


0.6143 ± 0.0050 


3.0 


0.003 


0.7544 ± 0.0005 


0.0775 ± 0.0013 


0.0815 ± 0.0015 


3.0 


0.02 


0.4724 ± 0.0040 


0.5162 ± 0.0021 


0.6602 ± 0.0053 


4.0 


0.01 


0.8460 ± 0.0002 


0.1459 ± 0.0008 


0.1631 ± 0.0010 


4.0 


0.03 


0.6996 ± 0.0008 


0.4600 ± 0.0013 


0.5554 ± 0.0016 



TABLE V: Real and imaginary part of the estimated fixed 
point J for the complex polarization, and corresponding dipole 
moment jj, per H2 unit, for periodic hydrogen chains with 
bond- length alternation L and applied electric field £. 



L(a. u. 


) DMC 


MP3(a) MP3(b) MP4(a) 


MP4(b) 


2.5 


50.57 ± 0.50 


51.35 54.33 50.02 


53.56 


3.0 


27.03 ± 0.55 


25.66 27.01 24.94 


26.51 


4.0 


16.04 ± 0.10 


15.4 16.13 15.00 


15.83 



TABLE VL DMC linear susceptibility a in a. u. per H2 unit 
for periodic linear hydrogen chains with bond-length alterna- 
tion L, compared with quantum chemistry results from Ref. 
[3| reported for the basis set (6)-31G(*)* (a) and (6)-311G(*)* 
(b). 



L(a. u. 


) DMC 


MP3(a) MP3(b) MP4(a) MP4(b) 


2.5 


652. ± 61. 




3.0 


89.8 ± 12. 


65.73 ± 0.06 73.00 ± 0.05 65.77 ± 0.11 74.68 ± 0.05 


4.0 


16.5 ± 1.1 





TABLE VII: DMC second hyper-susceptibility 7 in units of 
10'^ a. u. per H2 unit for periodic linear hydrogen chains 
with bond-length alternation L; for the chain with L—3.0 
a. u. , quantum chemistry MP results are taken from Ref. 
01, and are reported for the basis sets (6)-31G(*)* (a) and 
(6)-311G(*)* (b). 





Difference 


Relative Difference 


L(a. u. ) 


1 a 7 


a 


7 


2.5 


-3.69 260.70 


-7.3 


40. 


3.0 


-1.47 33.8 


-5.4 


37. 


4.0 


-1.35 3.63 


-8.4 


22. 



TABLE VIII: Differences between DMC and HP values for the 
linear susceptibility a in a. u. per H2 unit and for the second 
hyper-susceptibility 7 in 10"^ a. u. for H2 unit and relative 
differences in percent, for hydrogen chains with bond-length 
alternations L. 
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I/(a. u. 


) £ (a. u. ; 


) (a. u. ) 


2.5 


0.003 


2.396 ± 0.009 


2.5 


0.01 


2.873 ± 0.035 


3.0 


0.003 


1.751 ± 0.005 


3.0 


0.01 


2.261 ± 0.040 


4.0 


0.01 


1.391 ± 0.003 


4.0 


0.03 


1.620 ± 0.014 



TABLE IX: DMC localization along the longitudinal di- 
rection of hydrogen chains with bond-length alternations L 
for electric fields £. 



N 


a 


7 


10 


16.04 ± 0.10 


16.5 ± 1.1 


20 


15.70 ± 0.17 


12.7 ± 2.9 



TABLE X: DMC linear susceptibility a in a. u. and second 
hyper-susceptibility 7 in units of 10^ a. u. per H2 unit for 

periodic linear hydrogen chains with bond-length alternation 
L=4.0 a. u. calculated with supercells containing A'' H2 units. 



S (a. u. ) jRjzDMc} 5{^DMc} Mdmc (a. u. ) ^{z^y} "^{znY} Mhf (a. u. ) 

0. 0.6146 ± 0.0011 0.1249 ± 0.0019 0.1436 ± 0.0019 0.6311 0.0000 0. 

0.003 0.6127 ± 0.0008 0.1341 ± 0.0013 0.1544 ± 0.0018 0.6118 0.1443 0.1654 

0.01 0.6102 ± 0.0011 0.1608 ± 0.0024 0.1846 ± 0.0029 0.4021 0.4536 0.6055 

TABLE XI: Longitudinal dipole moment n and fixed point a calculated through DMC and HF for an hydrogen chain with 
bond-length alternation of 2.5 a. u. and for an electric field of 0.003 a. u. The trial wavefunctions are obtained from HF 
calculations with electric fields £. 



